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Abstract 

The basic optimal transportation problem consists in finding the 
most effective way of moving masses from one location to another, 
while minimizing the transportation cost. Such concept has been 
found to be useful to understand various mathematical, economical, 
and control theory phenomena, such as Witsenhausen's counterexam- 
ple in stochastic control theory, principal-agent problem in microeco- 
nomic theory, location and planning problems, etc. 
In this work, we focus on mobile association problems: the determina- 
tion of the cells corresponding to each base station, i.e., the locations 
at which intelligent mobile terminals prefer to connect to a given base 
station rather than to others. This work combines game theory and 
optimal transport theory to characterize the solution based on fluid 
approximations. We characterize the optimal solution from both the 
global network and the mobile user points of view. 
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1 Introduction 



Future wireless networks will be composed by intelligent mobile terminals, 
capable of accessing multiple radio access technologies, and able to decide for 
themselves the wireless access technology to use and the access point to which 
to connect. Within this context, we study the mobile association problem, 
where we determine the locations at which intelligent mobile terminals prefer 
to connect to a given base station rather than to others. We consider that 
these capabilities should be taken into account in the design and strategic 
planning of wireless networks. We analyze the case where mobile terminals 
within a cell share the same spectrum, and consequently, mobile terminals' 
decisions to which base station to connect affects the decision process of the 
other mobile terminals in the network. From these interactions, mobile ter- 
minals learn their optimal access point (self- learning), where the optimality 
depends upon the context. 

Starting from the seminal paper of Hotelling [T] a large area of research 
on location games has been developed. In PQ, the author introduced the 
notion of spatial competition in a duopoly situation. Plastria [2] presented 
an overview of the research on locating one or more new facilities in an envi- 
ronment where competing facilities already exist. Gabszewicz and Thisse j3] 
provided another general survey on location games. Altman et al. [I] studied 
the duopoly situation in the uplink scenario of a cellular network where the 
users are placed on a line segment. The authors realized that, considering 
the particular cost structure that arises in the cellular context, complex cell 
shapes are obtained at the equilibrium. Our work focuses on the downlink 
scenario and in a more general situation where a finite number of base sta- 
tions can compete in a one-dimensional and two-dimensional case without 
making any assumption on the symmetry of the users location. In order to 
do that, we propose a new framework for the mobile association problem 
using optimal transport theory (see [5] and references therein). This the- 
ory was pioneered by Monge [6] and Kantorovich [7\ and it has been proven 
to be useful in many mathematical, economical, and control theory con- 
texts [HI El [TD] - There is a number of works on "optimal transport" (see [TT] , 
and references therein) however the authors in [11] consider only an optimal 
selection of routes but do not use the rich theory of optimal transport. The 
works on stochastic geometry are similar to our analysis of wireless networks 
(see e.g. [T2] and references therein) but in our case we do not consider any 
particular deployment distribution function. Fluid models allow us to have 
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this general deployment distribution function. 

In the current work, we determine the spatial locations at which intelli- 
gent mobile terminals would prefer to connect to a given base station rather 
than to other base stations in the network. We obtain as well the spatial loca- 
tions which are more convenient from a global and centralized point of view. 
Obviously, in both approaches the optimality depends upon the context. In 
both considered cases, our aim is the minimization of the total power of 
the network, which can be considered as an energy-efficient objective, while 
maintaining a certain level of throughput for each user connected to the net- 
work. We propose the rich theory of optimal transport as the main tool of 
modelization of these mobile association problems. We are able to charac- 
terize these mobile associations under different policies and give illustrative 
examples of this technique. 

The remaining of this paper is organized as follows. Section [2] outlines 
the problem formulation of minimizing the total network power under quality 
of service constraints. We address the problem for the downlink case. Two 
different policies are studied: round robin scheduling policy (also known as 
time fair allocation policy) and rate fair allocation policy defined in Sec- 
tion [2] and studied in detail in Section 0] and Section [6] with uniform and 
non- homogeneous distribution of users. In Section [8] we give numerical exam- 
ples for one-dimensional and two-dimensional mobile terminals deployment 
distribution functions. Section M concludes the paper. 



2 The system model and problem formula- 
tion 

A summary of the notation used on this work can be found in Table [TJ 
We consider a network deployed on a region, denoted by T>, over the two- 
dimensional plane. The mobile terminals (MTs) are distributed according to 
a given deployment distribution function f(x,y). To fix ideas, if the consid- 
ered region is a square T> = [0 km, 1 km] x [0 km, 1 km] and the distribution 
of the users is uniform f(x,y) = 1, then the proportion of users in the sub- 
region A = [0 km, 1/2 km] x [0 km, 1/2 km] is 

II A f( x ^y) dxd y _ Aiea(A) _ 1 
jj v f(x,y)dxdy Area(£>) 4' 
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Table 1: Notation 



N 


Total number of MTs in the network 


K 


Total number of BSs 


f 


Deployment distribution of MTs 




Position of the i-th BS 


Q 


Cell determined by the i-th BS 




Number of MTs associated to the i-th BS 


Mi 


Number of carriers offered by the i-th BS 




Penalization function of non-service 


hi 


Channel gain function in the i-th cell 


6 


Path loss exponent in the i-th cell 



4 



The first equality is obtained because the distribution of the users is 
uniform. However, the expression at the left-hand side is general and it is 
always equal to the proportion of mobiles in a sub-region A. To simplify the 
notation, we normalize the function / such that 

J J f(x,y)dxdy = 1. (2) 

Consequently, the function / is a measure of the proportion of users over the 
network. 

The number of MTs in a subset of the network area A C Z>, denoted 
by N(A), is given by 

N{A) = N (JJ f{x,y)dxd y y (3) 

where N is the total number of MTs in the network. The integral on the 
right hand side between brackets takes into account the proportion of MTs 
distributed over the network area A. 

Examples of distribution of the users f(x,y): 

1. If the users are distributed uniformly over the network, then the mea- 
sure of the proportion of the users is given by 

/(*,y) = |!p (4) 

where \T>\ is the total area of the network, and K is a coefficient of 
normalization so that equation fl2]) holds. In this particular case, K = 1. 

2. If the users are distributed according to different levels of population 
density, then the measure of the proportion of the users would be 

{/hd if {%, y) is at a High Density region, 
/nd if (x, y) is at a Normal Density region, (5) 
/ld if (%■> y) is at a Low Density region, 

where /hd > /nd > /ld are defined similarly to equation (j3J with 
constants of normalization K w , K ND , K LD , such that Kn D > K ND > 
KjjD- As a particular example, we could consider in the domain [0, f] x 
[0, f] as High Density region the area [0, 1/4) x [0, 1] with Khd = 2, as 
Normal Density region the area [1/4, 1/2) x [0, 1] with i^ND = 1, and 
as Low Density region the area [1/2, 1] x [0, 1] with K^n = 1/2. In that 
case, equation holds. 
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3. If the distribution of the users is radial with more mobile terminals in 
the center of the network area and less mobile terminals in the suburban 
areas then 

f(x,y) = t- , (6) 

where Rd is the radius of the network and K D is a coefficient of nor- 
malization. 

Notice that the distribution of users f(x,y) considered in our work is more 
general than all the examples mentioned above. 

In the network, we consider K base stations (BSs), denoted by BSj ,i G {1, 
located at the fixed positions (xi,yi) , i e {1, . . . ,K}. The interference be- 
tween the different BS signals is ignored. We assume that the neighbouring 
BSs transmit their signals in orthogonal frequency bands. Furthermore, we 
assume that interference between BSs that are far from each other is negli- 
gible. Consequently, instead of considering the SINR (Signal to Interference 
plus Noise Ratio), we consider as performance measure the SNR (Signal to 
Noise Ratio). We consider the downlink case (transmission from base sta- 
tions to mobile terminals) and assume that each BS is going to transmit only 
to MTs associated to it. We denote by Cj the set of mobiles associated to 
the i-th BS, and by iVj the number of mobiles within that cell, both quan- 
tities to be determined. Notice that since the distribution of users f(x,y) 
considered in our work is general, instead of considering a particular distri- 
bution of mobiles, that we denote f(x,y), and an average throughput, that 
we denote 9(x,y), in each location (x,y), we can consider a constant average 
throughput 9 > and we can vary the distribution of mobiles f(x,y) such 
that the following equation holds: 

f(x,y)9 = f(x,y)9(x,y). (7) 

This would simply translate in the fact that for example mobile terminals 
with double demand than others would be considered as two users with the 
same demand. This can be done because of the fluid approximation of the 
network. 

If the number of mobiles is greater than the maximum number of carriers 
available in the i-th cell, denoted by MAXj, we consider a penalization cost 
function given by 

rAn _J XNi <MAXi, (R 

Ki[i) \ RiiNi-MAXi) if TVj > MAXj. { ) 
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We assume that be either a constant or a non- decreasing function^. 

We first study the case iVj < M and we study the general case in Section 16.11 
The power transmitted from BSj to an MT located at position (x, y), is de- 
noted by Pi(x, y). The received power at an MT served by BSj is Pi(x, y)hi(x, y). 
We shall further assume that the channel gain corresponds to the path loss 
given by hi(x, y) = (\/R 2 + d 2 (x, y))~^ where £ is the path loss exponent, R is 
the height of the base station, and di(x, y) is the distance between a MT at po- 
sition (x, y) and BSj located at (xi,yi), i.e., di(x,y) = ■sj (xi — x) 2 + (yt — y) 2 . 
The SNR received at mobile terminals at position (x, y) in cell Ci is given 
by SNRj(x, y) = Pi(x,y)hi(x,y)/a 2 , where a 2 is the noise power. We as- 
sume that the instantaneous mobile throughput is given by the following 
expression, which is based on Shannon's capacity theorem: 

Ofay) = log(l + SNR i {x,y)). (9) 

We want to satisfy an average throughput for MTs located at posi- 
tion (x,y) given by 9(x,y) > 0. We shall consider for this objective two 
policies defined in [13J: 



(A) Round robin scheduling policy: where each BS devotes an equal frac- 
tion of time for the transmission to each MT associated to it, and 

(B) Rate fair allocation policy: where each base station BS maintains a 
constant power sent to the mobile terminals within its cell and modifies 
the fraction of time allowed to mobile terminals with different channel 
gains, such that the average transmission rate demand is satisfied. 

For more information about this type of policies in the one dimensional 

case, see 



2.1 Round robin scheduling policy: Global Optimiza- 
tion 

Following this policy, BSj devotes an equal fraction of time for transmission 
to MTs located within its cell. The number of MTs located in the z-th cell 
is iVj, to be determined together with the cell boundaries. As BSj divides its 
time of service proportional to the quantity of users within its cell, then the 

1 For example, the maximum number of possible carriers in WiMAX is around 2048, so 
by using this technology we have MAX; = 2048. 
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throughput is given by 9f R (x, y) = 9i(x, y)/Ni. In order to satisfy a through- 
put 9(x,y), 6^ R (x,y) > 9(x,y), or equivalently, in terms of the power 

Pi(x,y)>-^-(2 N ^-l). (10) 

As our objective function is to minimize the total power of the network, the 
constraint will be reached, and we obtain by replacement 

Pfay) = a 2 (2 N ^ - 1)([R 2 + d^,y)] 1/2 ) +e - (11) 

From last equation we observe that: a) if the quantity of mobile terminals 
increases within the cell, the base station will need to transmit more power 
to each of the mobile terminals. The reason to do that is because the base 
station is dividing each time-slot into mini-slots with respect to the number of 
the mobiles within its cell, and b) the function (R 2 + d 2 (x, y))^ 2 on the right 
hand side give us the dependence of the power with respect to the distance 
between the base station and the mobile terminal located at position (x,y). 

Our objective is to find the optimal mobile association in order to mini- 
mize the total power of the network. Then as the total power 

K 

Ptota^^n^, (12) 

where P™** = \\ P{x, y)f(x, y) dx dy, (13) 

JJCi 

pmtra j g ^ e intracell power consumption in cell Cj. 

The global optimization for the mobile association problem is to deter- 
mine the cells Cj, % G {1, . . . , K}, to minimize the total power of the network: 

(RR) minV // P l (x,y)f(x,y)dxdy, (14) 

i=l JJ Ci 

subject to (ITT]), where f(x,y) is the deployment distribution function of the 
users. We solve this problem in Section HJ 

2.2 Generalized a-fairness formulation for minimiza- 
tion problems 

The general formulation for the problem of maximization of a function of 
the throughput given the constraint on the maximal power used admits a 
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generalized a-fairness formulation given by: 

max Ei^W"- 1 ]- ( 15 ) 

where we can identify different problems for different values of a: 

• a — maximization of throughput problem 

• a — > 1 proportional fairness (a uniform case of Nash bargaining) 

• a = 2 delay minimization 

• a — > +oo maxmin fairness (maximize the minimum throughput that a 
user can have). 

Since in minimization problems, the formulation is different since we are 
minimizing the total power on the network given the constraint of a mini- 
mum level of throughput, we define the following formulation, that we call 
generalized a-fairness for minimization problems: 

min E / — ^r[/( p ^)) a_1 - ?/) dx d y> (!6) 

where we can also identify different problems for different values of a: 

• a — maximization of the inverse of power (energy efficiency maxi- 
mization) 

• a — > 1 proportional fairness 

• ct = 2 minimization of total power 

• a — > +oo minmax fairness (to minimize the maximum power per BS). 

Note: the minmax fairness is not well studied in the literature but one 
can map the maxmin fairness studies into the minmax fairness for minimiza- 
tion problem. The convexity properties required becomes concavity, Schur 
convexity, sub-stochastic ordering, etc. 
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2.3 Rate fair allocation policy: User Optimization 

In the rate fair allocation policy, each BS will maintain a constant power 
sent to MTs within its cell, i.e., Pi(x, y) = Pi for each MT at location (x, y) 
inside cell Cj. However, the BS modifies the fraction of time allotted to 
MTs, set in such a way that the average transmission rate to each MT with 
different channel gain is the same, denoted by 6(x, y), for each mobile located 
at position (x, y). 

We study the equilibrium states where each MT chooses the BS which will 
serve it. Given the interactions with the other mobile terminals it doesn't 
have any incentive to unilaterally change its strategy. A similar notion of 
equilibrium has been studied in the context of large number of small players 
in road-traffic theory by Wardrop [H] . 

Definition.- The Wardrop equilibrium is given in the context of cellular 
systems by: 



A Wardrop equilibrium is the analog of a Nash equilibrium in the case of 
a large number of small players, where, in our case, we consider the mobile 
terminals as the players. In this setting, the Wardrop equilibrium indicates 
that if there is a positive proportion of mobile terminals associated to the 
i-th base station (the left-hand side condition in (I17al) ). then the throughput 
that the mobile terminals obtain is the maximum that they would obtain 
from any other base station (right-hand side consequence in (I17al) ). The 
second condition indicates that if there is one base station that doesn't have 
any mobile terminal associated to it (left-hand side condition in fll7al) ). it is 
because the mobile terminals can obtain a higher throughput by connecting 
to one of the other base stations (right-hand side consequence in (1 1 7a.j) ) . 

We assume that each base station is serving at least one mobile terminal, 
(if that is not the case, we remove the base station that is not serving any 
mobile terminal). Then, the equilibrium situation is given by 



To understand this equilibrium situation, consider as an example the sim- 
ple case of two base stations: BSi and BS2. Assume that BSi offers more 




(17a) 




(17b) 



6\ = 6 2 



(18) 



10 



throughput than BS2. Then, the mobile terminals being served by BS2 will 
have an incentive to connect to BSi. The transmitted throughput depends 
inversely on the quantity of mobiles connected to the base station. As more 
mobile terminals try to connect to base station BSi the throughput will di- 
minish until arrive to the equilibrium where both base stations will offer the 
same throughput. 

The condition 9\ = 9 2 = ■ ■ ■ = 9k is equivalent in our setting to the 
condition t\ = r 2 = . . . = rx, and also to the condition SNRi = SNR 2 = 
. . . = SNRr-. Let us denote by r to the rate offered by the base station at 
equilibrium, i.e., r := r\ = r 2 = ... = tk and we denote by (3 to the SNR 
offered by the base station at equilibrium, i.e., (3 = SNR4 = SNR 2 = . . . = 
SNR^. 

This condition in terms of the power is equivalent to 

P^ = p[[ a 2 ([R 2 + d 2 (x,y)} 1 / 2 )^f(x,y)dxdy. (19) 

JJCi 

We want to choose the optimal mobile assignment in order to minimize 
the total power of the network under the constraint that the mobile terminals 
have an average throughput of 9, i.e., 

K 

MinP total = ^^ intra . (20) 
i=i 

Then our problem reads 

(RF) Min^ [[ a 2 (R 2 + d 2 (x,y)f 2 f(x,y)dxdy. (21) 
c * <=i JJ c> 

We will solve this problem in Section El Thanks to optimal transport 
theory we are able to characterize the partitions considering a general setting. 
In the following section, we will briefly describe optimal transport theory 
and motivate the solution of the previously considered mobile association 
problems. 

3 Basics in Optimal Transport Theory 

The theory of mass transportation, also called optimal transport theory, goes 
back to the original works by Monge in 1781 [6], and later in 1942 by Kan- 
torovich [7]. The work of Brenier [?] has renewed the interest for the subject 
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and since then many works have been done in this topic (see [5] and references 
therein) . 

The original problem of Monge can be interpreted as the question: 

"How do you best move given piles of sand to fill up given holes 
of the same total volume?" . 

In our setting, this problem is of main importance. Suppose that base stations 
are sending information to mobile terminals in a grid area network and the 
positions of base stations and mobile terminals are given. 

What is the "best move" of information from the MTs to the 
BSs? 

Both questions share similarities as we will see. The general mathematical 
framework to deal with this problem is a little technical but we encourage to 
focus on the main ideas. 

The framework is the following: 

We first consider a grid area network T> in the one- dimensional case. As 
an example, the function f(t) will represent the proportion of information 
sent by mobile terminals 

d(i(t) := f(t) dt. (22) 

The function g(s) will represent the proportion of information received by a 
base station at location s 

du{s) := g(s)ds. (23) 

The function T (called transport map) is the function that transfers infor- 
mation from location s to location t. It associates mobile terminals to base 
stations and transports information from base stations to mobile terminals. 
Then the conditions that each mobile terminal satisfies its downlink demand 
is written 

/ F(y)g(y)dy= [ F(T(x))f(x)dx, (24) 

J A J{x£X:T(x)£A} 

for all continuous function F, where X is the support^] of function / and we 
denote this condition (following the optimal transport theory notation) as 

T#/i = v, (25) 

2 The support of a function / is the closure of the set of points where the function is 
not zero, i.e., support (/) = {t : f(t) ^ 0} 
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which is an equation of conservation of the information. Notice that, in 
communication systems there exists packet loss so in general this constraint 
may not be satisfied, but considering an estimation of the packet loss by 
sending standard packets test, this constraint can be modified in the reception 
measure v. If we can not obtain a good estimation of this reception measure, 
we can consider it in its current form as a conservative policy 

In the original problem, Monge considered that the cost of moving a 
commodity from position x to a position y depends on the distance c(\x — y\). 
Then the cost of moving a commodity from position x through T to its new 
position T{x) will be c(\x — T(x)\). For the global optimization problem, we 
consider the additive total cost over the network, which in the continuum 
setting will be given by 

Min / c(\x — T(x)\) f(x) dx such that T#fi = v, (26) 
Jv 

where /i and v are probability measures and T : T> — > T> is an integrable 
function. This problem is known as Monge's problem in optimal transport 
theory. 

The main difficulty in solving Monge's problem is the highly non-linear 
structure of the objective function. For examples on the limitations on 
Monge's modelization, see [5]. As an example, consider the domain T> = 
[0,2], the transmission from a base station located at position 1, denoted 
fi — 8\, and the throughput demanded to this base station by two mobile 
terminals located at positions and 2, denoted v = |<5 + ^<5 2 - According to 
the formulation given by Monge, there is no splitting of throughput, since 
everything that is transmitted from one location has to go to another lo- 
cation. So this simple problem doesn't have a transport map (see Fig. [TJ). 
This limitation is due in part to the original considered problem, but as we 
will see this limitation is overcome by Kantorovich's approach. We have also 
pointed out the limitations of Monge's problem that motivated Kantorovich 
to consider another modeling of this problem in [TJ. 

Kantorovich noticed that the problem of transportation from one location 
to another can be seen as "graphs of functions" (called transport plans) in 
the product space (See Fig. [2]). 

The idea is to minimize the objective function over the space of graphs 
g = (Id x T) in the product space. Then with the condition that each mobile 
terminal satisfies its uplink demand and that the information is received at 
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MTl 



BS 



MT2 



Figure 1: Monge's problem can not model a simple scenario of two mobile 
terminals and one base station. Kantorovich's problem however can model 
very general scenarios. 

the base stations, Kantorovich's problem reads 



is denoted the ensemble of transport plans g, tti (x, y) stands for the projection 
on the first axis x, and 7r 2 (x, y) stands for the projection on the second axis y. 

The relationship between Monge and Kantorovich problems is that every 
transport map T of Monge's problem determines a transport plan g? = 
(Id x T)#/i in Kantorovich's problem with the same cost (where Id denotes 
the identity). However, Kantorovich's problem considers more functions than 
the ones coming from Monge's problem (which can always be viewed as the 
product of the identity and the map T), so we can choose from a bigger set 



Then, every solution of Kantorovich's problem is a lower bound to Monge's 
problem, i.e., 




(27) 



where 



n O> v ) = {9 ■ 7ri#9 = At and n 2 #g = v}, 



(28) 



n(/i, v). 



Min 



c(x,y)dg(x,y) 



< Min 

T#»=u 



[ c(\x-T(x)\)f(x)dx. (29) 



VxV 



JV 
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Figure 2: Kantorovich considered "graphs" where the projection in the first 
axis coincide with the mobile terminal position (MTi = 3.5, MT2 = 5 and 
MT 3 = 6.5) and the second axis coincides with the base station position 
(BSi = 4 and BS 2 = 6). 



Theorem 3.1 Consider the cost function c(\x—y\) = \x— y\ p . Let fi andv be 
probability measures in T> and fix p > 1. We assume that /i can be written^ as 
dfj, = f(x) dx. Then the optimal value of Monge's problem coincides with the 
optimal value of Kantorovich 's problem, i.e., M p (fjL,v) = W p (/i, v) and there 
exists an optimal transport map from fi to v, which is also unique almost 
everywhere if p > 1. 

This result is very difficult to obtain and it has been proved only recently 
(see [?] for the case p — 2, and the references at [5] for the other cases). 

The case that we are interested in can be characterized because the image 
of the transport plan is a discrete finite set. 

Thanks to optimal transport theory we are able to characterize the par- 
titions considering general settings. To this purpose, consider locations 
(xi,yi) . . . , (xk, yic)i the Euclidean distance di(x, y) = yjjx — Xi) 2 + (y — yi) 2 , 
and F a continuous function. 

3 The exact condition is that /i is absolutely continuous with respect to the Lebesgue 
measure. A probability measure n is absolutely continuous with respect to the Lebesgue 
measure if the function F{x) = 00, x\) is locally an absolutely continuous real function. 
A function / is an absolutely continuous real function if there exists an integrable function g 
such that /0) = /(a) + g(t) dt 
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Theorem 3.2 Consider the problem 

K 



(PI) mnf^jj^ F(d l (x,y)) + s l ^JJ^f(io,z)doodz 



f(x,y) dx dy, 
(30) 



where Ci is the cell partition ofT>. Suppose that are continuously differ- 
entiable, non- decreasing, and convex functions. The problem (PI) admits a 
solution that verifies 



(31) 



C t = {(x,y):F(d i (x,y)) + s i (N i ) + N i -s' i (N i 
(SI) { < F{d 3 (x,y)) + s 3 {N 3 ) + Nj ■ s'^N,)} 

Ni = JJ C f(co,z) du dz. 

Proof.- See Appendix A 
Theorem 3.3 Consider the problem 



(P2) MmJ2 JJ c F(d z (x,y))-m l ^JJ^f(u,z)dudz^ f(x,y)dxdy, 

(32) 

where Ci is the cell partition of V. Suppose that rrii are derivable. The 
problem (P2) admits a solution that verifies 



(S2) 



Ci = {(x,y) : mi(Ni)F(di(x,y)) f(x,y) + Ui(x,y) 

<m j (N j )F(d j (x,y))f(x,y) + U j (x,y)} 
Ui = m'^Ni) ff c . F(di(x, y))f{x, y) dx dy 
Ni = JJ C f(u, z) du dz. 



(33) 



Proof.- See Appendix B 

Notice that in problem (PI) if the functions = the solution of the 
system (SI) becomes the well known Voronoi cells. In problem (P2) if we 
have that the functions hi = 1 we find again the Voronoi cells. In general 
however, the Voronoi configuration is far from being optimal. 
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Interference vs location when BS1 is at and BS2 is at -10 




Location of mobile terminals 



Figure 3: Interference as a function of location of mobile terminals when BSi 
is at position (solid line) and BS 2 at —10 (dotted line). 



SINR vs location when BS1 is at and BS2 is at -10 




Location of mobile terminals 



Figure 4: SINR as a function of location of mobile terminals when BSi is at 
position (solid line) and BS 2 at —10 (dotted line). 
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Figure 5: Wardrop Equilibrium: Thresholds determining the cell boundaries 
as a function of the location of the base stations. The network is deployed over 
the interval [—10, 10] (one-dimensional case) presented here vertically. We 
consider a uniform distribution of MTs and we find the threshold (solid line) 
determining the cell boundaries as a function of the base stations positions 
(dotted lines) by changing the position of one of them. BSi is fixed at 
position and we change the position of BS2 from —10 to +10. 
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Figure 6: Worst Equilibrium: Thresholds determining the cell boundaries 
(vertical axis) that give the worst equilibrium in terms of the SINR as a 
function of the location of BS 2 for BSi at position 0. 



4 Round Robin Scheduling Policy 

We assume that a service provider wants to minimize the total power of 
the network while maintaining a certain average throughput of 9 to each 
mobile terminal of the system using the round robin scheduling policy given 
by problem 

(RR) MinV ff a 2 (R 2 + d i (x,y) 2 f /2 (2 N > 9 -l)f(x,y)dxdy. (34) 

We see that this problem is an optimal transportation problem like the 
one in (PI) with cost function given by 

F(d l (x,y)) = ( r 2 (R 2 + d t (x,yff 2 and (35) 
mi(x,y) = (2 N *-1). (36) 

From the previous theorem, we can derive an explicit expression for this 
configuration. 
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Proposition.- There exist a unique optimum given by 
Q ={x e V : di(x , yo) p + hiiNi) + A^(^) 
< d^yoY + kj(Nj) + NtfiNj), Vj ^ 1} 

Ni = J J f(x , y ) dx dy . 
Let's see a direct application of our results: 



(37) 
(38) 



Example 4.1 Consider a network of N = 2500 mobile terminals distributed 
according to f(x) in [0, L] (for example, with L = 5.6 miles for WiMAX 
radius cell). We consider two base stations at position BSi = and BS2 = L 
and R—l. Then, the mobile association threshold (the boundary between 
both cells, i.e., the location at which the mobile terminals obtain the same 
throughput by connecting to any of both base stations) is reduced to find x 
such that the following equality holds: 



(2 



N 2 e 



{2 Nlt) - l)(l + x 2 )f{x) + 2™ 6 log 2 
l)(l + (l-x) 2 )/(x) + 2 JV2 ^log2 



x 



4 9 X 6 

--2x + x'- — 



(39) 



Notice that this is a fixed point equation on x. If the mobile terminals 
are distributed uniformly, the optimal solution is given by C\ = [0,1/2 L) 
and C2 = [1/2 L,L], which is the solution that Voronoi cells would give 
us and in that case the number of mobile terminals connected to each base 
station would be given by 



Nx = N 2 = 1250. 



(40) 



However, if the deployment distribution of the mobile terminals is more con- 
centrated near BS 2 than BSi, consider for example f(x) = 2x, the optimal 
solution is given by C\ = [0, q) and C 2 = [q, L] with q = 0.6027 L and 



^ = 908 and N 2 = 1592. 



(41) 



Notice that in the global optimization solution, the number of mobile terminals 
connected to BSi is smaller that the number of mobile terminals connected 
to BS2. However, the cell size is bigger. 
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(a) Equilibrium in Non-Uniform distribution 

Figure 7: Thresholds determining the cell boundaries (vertical axis) of the 
best equilibrium in terms of the SINR as a function of the location of BS2 for 
BS1 at position when we consider a non- homogeneous distribution given by 
f(x) = (L — x)/2L 2 . Example: equilibrium when then distribution of mobile 
terminals is given by f(x) = 2x in the interval [0, L] and the positions of the 
base stations are BSi = and BS2 = L. 
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Figure 8: Wardrop Equilibrium in the Non- homogeneous case: The network 
is deployed over the interval [—10, 10] (one-dimensional case) presented here 
vertically. We find the threshold (solid line) determining the cell boundaries 
as a function of the base stations positions (dotted lines) by changing the 
position of one of them. BSi is fixed at position and we change the position 
of BS2 from —10 to +10. The deployment distribution of the MTs is given 
by f(x) = (L-x)/2L 2 . 
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5 Fairness problem 



As we mention in section [2J the solution given by previous section H] is optimal 
but may not be fair to all the mobile terminals since it will give higher 
throughput to the mobile terminals that are near the base stations. To deal 
with this problem we considered the fairness problem given by 

Min E JJ c ^ 1 (v\R 2 + dKx,y))' / r~ 1 (2 Nie -V a - 1 f^y)dxdy. 

As we can see this is also an optimal transportation problem (PI) where the 
functions considered in this setting are given by 

F(di(x,y)) = -J- {^(R' + d^y))^ 1 (42) 
a — 1 

m t (x,y) = (2 N > e -ir~ 1 . (43) 

Using Theorem 13.21 we are able to characterize the optimal cells for any a 
considered. 



6 Rate fair allocation policy 

In this framework, we give the possibility to mobile terminals to associate to 
the base station they prefer in order to minimize their power cost function 
while maintaining, as quality of service measurement, an average throughput 
of 6. 

As we presented in Section [21 this problem is equivalent to 

(RF) MinV // a 2 (R 2 + d 2 (x,y)f 2 f(x,y)dxdy. (44) 
c ' »=i JJCi 

Notice that this problem is equivalent to (PI) where the functions 

F(d t (x,y))=(a 2 (R 2 + d t (x,y) 2 ^ 2 ), (45) 
and Sj = 1. The problem has then a solution given by 
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Proposition.- There exist a unique optimum given by 

Q= {xeV:a 2 (R 2 + d*(x ,y )y/ 2 

<a 2 (R 2 + d 2 (x ,y ))Z/ 2 ,Vj^t} (46) 

N i = IIcJ^Vo) dx d V0 > 

which is represented by the Voronoi cells. 



6.1 Penalization function 

Notice that the penalization function or the case when the number of mobile 
terminals is greater than the number of carriers available in the cell in the 
rate fair allocation policy case is equivalent to (PI) where the functions 

F(di(x, y)) = {cr 2 (R 2 + d^x, yff' 2 ) . (47) 



(N\- (M\-S if^<MAX i? 

s^i) - Hi{JSi) - | _ MAX ^ . f N _ > MAX ^ (48) 

The problem has then a solution given by 
Proposition.- There exist a unique optimum given by 

d = {(x,y) : a 2 (R 2 + d 2 (x ,y ))^ 2 + s^N,) + N, ■ s'^N,) 

< a 2 (R 2 + d 2 (x ,y ))V 2 + s^Nj) + N, ■ s'^N,)} (49) 
Ni = ff c f(ou,z) duo dz. 



7 Performance Gap 

As a second example consider again the case when the mobile terminals are 
uniformly distributed on V = [0,1] but this time the two antennas are located 
at coordinates X\ — and x 2 — 1. Consider the case when p — 1 and 

n lfl(1 , , , J for 0<x< 0.999 
Sl (x) = 100 and *,(*) = | x for q.999 < £ < 1. 

Then the equilibrium cell configuration (Cf , Qf ) is given by 
Cf = and Cf = [0,l], 
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Cell Boundaries with uniform distribution 
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Figure 9: Wardrop Equilibrium with Multiple Base Stations: The grid area 
network is the interval [—10, 10] presented here vertically. We consider a 
uniform distribution of MTs. We find the threshold (solid and dashed lines) 
determining the cell boundaries as a function of the base stations positions 
(dotted lines) by changing the position of one of them. BSi is fixed at 
position —10 and BS2 is fixed at position 10 and we change the position of 
BS 3 from -10 to +10. 

and the optimum cell configuration (C*,^) is 

C* = [0,001[ and C 2 * =] 0.001, 1]. 

The optimum is very unfair for mobile terminals living in the first cell C*, 
who pay x + 100, whereas the other mobile terminals just pay the distance 
from 1. This is a toy example but it gives an idea of the performance gap 
between the centralized and the decentralized scenarios, also known as Price 
of Anarchy. 

8 Numerical Simulations 

In this section, we present several numerical results that validate our theo- 
retical model. 
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Figure 10: Wardrop Equilibrium in the 2D case: The grid area network is the 
square [—4,4] x [—4,4]. We set the noise parameter a = 0.3 and we set four 
BSs at positions BSi = (-3,-3) BS 2 = (3,-3) BS 3 = (3,3) BS 4 = (-3,3) 
and one at the origin BS 5 = (0,0). We determine the cell boundaries (deep 
lines) for the uniform distribution of users 




Figure 11: 2D case: Cell contours of the best equilibrium with uniform 
distribution of users. 
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Figure 12: Wardrop Equilibrium in 2D Non-Uniform Case: The grid area 
network is the square [—4,4] x [—4,4]. We set the noise parameter a = 0.3 
and we set four BSs at positions BSi = (—3, —3) BS2 = (3, —3) BS3 = (3, 3) 
BS4 = (—3,3) and one at the origin BS5 = (0,0). We determine the cell 
boundaries (deep lines) for the non-uniform distribution of users given by 
f(x, y) = (L 2 — (x 2 + y 2 ))/K where K is a normalization factor. The latter 
situation takes into account when mobile terminals are more concentrated in 
the center and less concentrated in suburban areas. 
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Figure 13: 2D Non- Uniform: Cell contours of the best equilibrium with non- 
uniform distribution of users. 

8.1 One-dimensional case 

We first consider the one- dimensional case and we consider a uniform distri- 
bution of users in the interval [—10, 10]. We set the noise parameter a = 0.3. 
In Fig. |5l we fix one base station BS2 at position and take as parameter 
the position of base station BSi. We consider as path loss exponent £ = 2. 
Red lines shows the positions of the BSs. We are able to determine the 
cell boundary (solid blue curve) from BSi and BS2 at different positions. In 
Fig. |H1 we fix two base stations BSi = —10 and BS 2 = 10 and we take as 
parameter the position of base station BS3. Red lines shows the positions of 
the BSs. We determine the cell boundary (solid blue curve) from BSi and 
BS3 and the cell boundary (dotted blue curve) from BS2 and BS3. 

8.2 Two-dimensional case: Uniform and Non-Uniform 
distribution of users 

We consider the two-dimensional case. We consider the square [—4, 4] x [—4, 4] 
and the noise parameter a = 0.3. We set five base stations at positions BSi = 
(-3,-3), BS 2 = (3,-3), BS 3 = (-3,3), BS 4 = (3,3), and BS 5 = (0,0). 
We determine the cell boundaries for the uniform distribution of MTs (see 
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Fig. ITTj) and we compare it to the cell boundaries for the non-uniform distri- 
bution of MTs given by f(x, y) = (L 2 — (x 2 + y 2 ))/ K where if is a normal- 
ization factor. The latter situation can be interpreted as the situation when 
mobile terminals are more concentrated in the center and less concentrated 
in suburban areas as in Paris, New York or London. We observe that the 
cell size of the base station BS5 at the center is smaller than the others at 
the suburban areas. This can be explained by the fact that as the density 
of users is more concentrated in the center the interference is greater in the 
center than in the suburban areas and then the SINR is smaller in the center. 
However the quantity of users is greater than in the suburban areas. 

9 Conclusions and Future Perspectives 

In the present work, we have studied the mobile association problem in the 
downlink scenario. The objective is to determine the spatial locations at 
which mobile terminals would prefer to connect to a given base station rather 
than to other base stations in the network if they were offered that possi- 
bility (denoted decentralized scenario). We are also interested in the spatial 
locations which are more convenient from a centralized or from a network 
operator point of view. In both approaches, the optimality depends upon 
the context. In the considered cases, we consider the minimization of the 
total power of the network, which can be considered as an energy-efficient 
objective, while maintaining a certain level of throughput for each user con- 
nected to the network. We have proposed a new approach using optimal 
transport theory for this mobile association problems and we have been able 
to characterize these mobile associations under different policies. 

The present work can be extended in several different directions. One 
of these possible directions is to study the price of anarchy between the 
centralized and decentralized scenario. As we presented in Section [7J the 
considered example give us an indication that the price of anarchy should 
be unbounded but currently we don't have precise bounds. The price of 
anarchy should be studied in both scenarios: the sum of a function and 
the multiplication. It should be interesting to study the application in the 
particular case when the network is an LTE network. Since our model is quite 
simplified in order to obtain exact solutions we could include the cases for 
the fading and shadowing effects. It is implicitly considered that the number 
of users in the network is stationary, but since at different times of the day 



29 



there are different number of users, this management capabilities should be 
taken into account. 
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Appendix A 

Consider the problem (PI) 



where Cj is the cell partition of V. Suppose that are continuously differ- 
entiate, non-decreasing, and convex functions. The problem (PI) admits a 
solution that verifies 



Proof.- The proof is based on Proposition 3.5 of Crippa et al. [15]. We 
include the proof for completeness and because part of it (mainly the ex- 
istence of the solution) will be used in the proof of the following theorem. 
Notice that we have considered the case F(di(x, y)) = \(x, y) — (xj, yi)\ p , but 
this holds for any continuous function F. 

From Section [31 let us recall that Monge's problem can be stated as 
follows: given two probability measures, /i and u, and a constant p > 1 we 
consider the minimization problem (denoted by M p (yU, z/)): 





(51) 



M p (/i, v) := inf 




x — T(x)\ p dfi(x) 



i/p 



: T : T> — V T> Borel and such that T#/j, 



(52) 
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The relaxed formulation of Monge's problem (denoted by W p (/j,, v)) can 
be stated as follows 

W p (ji,u):= inf 1(1 \x-y\^dj(x,y)) }, (53) 
7en(At,v) \J VxV J J 

where IT(/i, i/) is the set of probability measures such that 771(7) = A 4 an d 
^2(7) — v where n\ is the projection on the first component and 7r 2 is the 
projection on the second component. 

If the probability measure fi can be written as dfi = f(x) dx (i.e. it 
is absolutely continuous with respect to the Lebesgue measure), then the 
optimal values of both problems coincide M p (/z, z/) = and there 

exists an optimal transport map from fi to v which is unique /-a.e. if p > 1. 
Another important characteristic of this relaxation is that it admits a dual 
formulation: 

\ f f , f , u(x) + v (y) < \x — y\ p for u-a.e. x and u-a.e. y 

(54) 

Moreover, there exists an optimal pair (u, v) for this dual formulation, and 
when u is an atomic probability measure (it can be written as u = biS yi ) 
the dual formulation becomes 

rir „, s If, \ , / x u(x) + v(yi) < \x — yA p for u-a.e. x and every i eN 

W^, u) = sup dp + g 6,„( W ) : " I e ™ v)jV £ Ll{v) 

(55) 

There exists another interesting characteristic when one of the measures is 
absolutely continuous with respect to the Lebesgue measure and the other 
measure is an atomic measure. If the probability measure \x can be written 
as d\i = f(x) dx where / is a nonnegative function, (yi)^ is a sequence of 
points in the domain V such that v = ^2 ie ^bi6 Vi , (Bi)^ is a partition of 
the domain such that the map T(x) = X^eN^lBi( x ) ^ s an optimal transport 
map from fi to u, the pair (u, v) is a solution of the dual formulation (JBl 
then 



u(x 



inf(|x - yi\ p - v(yi)) = }(\x - yi\ p - v{yi))l Bi {x) for f-a.e. x. (56) 



We also have a similar converse characteristic. If (Bi) ie ^ is a partition of X>, 
we set bi = f B . f(x) dx, v = Y,ieN b i^yi and T ( x ) = J2ieN Vi^Ax), and there 
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exists two functions u G L*(X>), v G L l u {V) satisfying the condition (15B"]) . 
then T is optimal for M p (fi,u) and the pair (w,t>) is optimal for the dual 
formulation fl54|) . 

In order to prove the existence and uniqueness of the solution we need to 
consider the following: We denote by S the unit simplex in M. k : 



S=<c=(cu...,c k )eR k . 
From Theorem I3.1[ we deduce that 

f(x) dx 



Ci > U, 2 

i=i 



(57) 



inf V 

partition of 7? 
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i=l 



(58) 

,c k ) = WP U, £)J =1 c^.J is 



The function F : S — > ffi. defined by -F(ci, . 

continuous and convex. 

If the functions Sj are lower semi-continuous, then there exists an optimum 
and if in addition the maps i]i(t) := thi(t) are strictly convex, the optimum 
is unique. 

We can then characterize the solution. If (Ai) i= i % is an optimum, and 
Si are differentiable in ]0, 1] and continuous in 0, then the following holds: 



Q ={xeV :\x- Xi\ p + Si(Ni) + NiS^Ni) 
<\x- xj\ p + s^Nj) + Njs'^Nj), Vj ^ i 



Ni 



f(x) dx. 



(59) 
(60) 



Appendix B 

Consider the problem (P2) 



F(di(x,y)) ■ rrii 



f(u, z) dui dz 



Ci 



f(x,y)dxdy, (61) 
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where Cj is the cell partition of V. Suppose that m 8 are derivable. The 
problem (P2) admits a solution that verifies 

r d= {x : m,(7v,)^R(^y))/(^y) + ^,y) 

, q9 x I < mjiNAFidfay)) f(x,y) + U 3 {x,y)} 

Ui= m>(N l )ff Ci F(d t (x,y))f(x,y)dxdy W 
y Ni — ff c f(u>, z) du> dz. 

Proof.- The proof is similar to Appendix A. The problem (P2) can be 
rewriten as follows: 

inf | W Ci dx, CiSx^J ■ > 0, °i = 1 J > ( 63 ) 

where 

W Ci fdx,y^Ci8 x . = inf ly^l F(\x -Xi\) ■m i (c i )f(x)dx : / / dx = q, (A) fe =1 partition 

(64) 

Let 7 be optimal for Monge-Kantorovich problem, then one can consider 

k 

l(x,y) = 7*0*0 <g> 8 x .(y), Ci = %(V), (65) 



where 7, is a positive measure such that < f dx 

The function F : S ^ R defined by F{c u ...,c k ) = (jj,, J2i =1 Q^) is 
lower semi-continuous. Then there exists a solution and it is unique almost 
surely. 

Let (cj)f =1 be a solution of (P2) and (v4«)f =1 the associated optimal par- 
tition. Let us fix two indices jo an d a point rr G A io . Let e > 0. Let us 
consider the open ball of radius e and center x that we denote as B t . We de- 
note its measure as c e . We make a small variation of the optimal partition by 
taking from A io the ball B e and adding it to Aj . Since the partition (A)f=i 
is optimal then 



\)h io (c io )f(x)dx + / F(| x Xj 1 ) hj (cj )f(x) dx 
< / F(\x - x io \)h io (c io - c e )f(x)dx + F(\x - x jo \)h jo (c jo + c £ )f(x) dx, 

JA in \B € JA in UB s 



(66) 
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which is equivalent to 

/ F(\x-x io \)(h io (c io )-hi (c io -c e ))f(x)dx+ / F(\x - x io \)h io (c io - c e )f(x) dx 
Ja H) JB t 

— / F{\x — Xi \)(h j0 {c jo + c e ) - h j0 (c j0 ))f(x) dx + / F(\x - x jo \)hj (cj + c £ )f(x) dx. 

JA jQ JB e 

(67) 

Dividing the previous equation by c e and taking the limit when e — > 0, 
we obtain 

F{\x — Xi \)K (ci )f(x) dx + F(\x - x io \)h io (c io )f(x ) 
< [ F(\x-x j0 \)h'Jc j0 )f(x)dx + F(\x-x j Xh j0 (c j0 )f(x ). (68) 

JA in 



Reorganizing the terms we obtain the desired result. 
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